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Abstract We report an improved measurement of the neutrino mixing angle #13 from the Daya Bay Reactor 
Neutrino Experiment. We exclude a zero value for sin 2 2#i3 with a significance of 7.7 standard deviations. 
Electron antineutrinos from six reactors of 2.9 GW t h were detected in six antineutrino detectors deployed 
in two near (flux-weighted baselines of 470 m and 576 m) and one far (1648 m) underground experimental 
halls. Using 139 days of data, 28909 (205308) electron antineutrino candidates were detected at the far hall 
(near halls). The ratio of the observed to the expected number of antineutrinos assuming no oscillations 
at the far hall is 0.944 ± 0.007(stat.) ± 0.003(syst.). An analysis of the relative rates in six detectors finds 
sin 2 26»i 3 = 0.089±0.010(stat.)±0.005(syst.) in a three-neutrino framework. 

Key words neutrino oscillation, neutrino mixing, reactor, Daya Bay 

PACS 14.60.Pq, 29.40.Mc, 28.50.Hw, 13.15.+g 
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1 Introduction 



until the analyses were finalized. 



Observations of neutrinos and antineutrinos pro- 
duced in the sun, the atmosphere, reactors, and 
from particle beams provide overwhelming evidence 
that the flavors of neutrinos change (oscillate) [ll-Q . 
The preponderance of data support a three-neutrino 
framework where three flavor states (v e ,iy^ 1 i> T ) are 
superpositions of three mass states (v 1 ,v 2 ,v:i)- This 
mixing can be quantified using a unitary 3x3 mix- 
ing matrix described in terms of three mixing angles 
(#12, #23, #13) and a CP violating phase (5) Neu- 
trino oscillations are also dependent on the differences 
in the squares of the neutrino masses. 

The Daya Bay collaboration recently measured a 
non-zero value for sin 2 20 13 = 0.092 ± 0.016(stat.) ± 
0.005(syst.) an observation consistent with previ- 
ous and subsequent experimental results 
absolute terms, the value of 9 13 is now known with 
better precision than either of the other two mixing 
angles. Constraining the value of f? 13 increases the 
constraints on the other mixing parameters (mixing 
angles and mass squared differences) through a global 
fit of all available oscillation data 12, 13^ . 

For reactor-based experiments, in a three- neutrino 
framework, an unambiguous determination of 9 13 can 
be extracted via the survival probability of the elec- 
tron antineutrino T> e at short distances (O(km)) from 
the reactors 

P sur « 1 - sin 2 26>i3 sin 2 ( 1.267 Am 2 , L/E) , (1) 

where Am 2 ! can be approximated by Am 2 tm = 



L3 
> L4 



AD3 Ling Ao-II NPP 



(2.32 



+0.12\ 



x 10" 3 eV 2 



14j , E is the V e energy in Me V 
and L is the distance in meters between the V e source 
and the detector (baseline). The near- far arrange- 
ment of antineutrino detectors (ADs), as illustrated in 
Fig. [TJ allows for a relative measurement by compar- 
ing the observed 77 e rates at various distances. With 
functionally identical ADs, the relative rate is inde- 
pendent of correlated uncertainties, and uncorrelated 
reactor uncertainties are minimized. 

The results reported here were derived using the 
same analysis techniques and event selection as our 
previous results [8| , but were based on data collected 
between December 24, 2011 and May 11, 2012, a 2.5 
fold increase in statistics. A blind analysis strategy 
was adopted for our previous results, with the base- 
lines, the thermal power histories of the cores, and 
the target masses of the ADs hidden until the analy- 
ses were finalized. Since the baselines and the target 
masses have been unveiled for the six ADs, we kept 
the thermal power histories hidden in this analysis 




' ' Daya Bay NPP 

Fig. 1. Layout of the Daya Bay experiment. 
The dots represent reactor cores, labeled as 
Dl, D2, LI, L2, L3 and L4. Six antineutrino 
detectors (ADs) were installed in three exper- 
imental halls (EHs). 

2 The Experiment 
2.1 Site 

The Daya Bay nuclear power complex is located 
on the southern coast of China, 55 km to the north- 
east of Hong Kong and 45 km to the east of Shen- 
zhen. A detailed descri ptio n of the Daya Bay exper- 
iment can be found in Tj^, As shown in Fig. 1, 
the nuclear complex consists of six reactors grouped 
into three pairs with each pair referred to as a nu- 
clear power plant (NPP). All six cores are function- 
ally identical pressurized water reactors, each with a 



maximum of 2.9 GW thermal power [ljj. The last 
core started commercial operation on Aug. 7, 2011. 
The distance between the cores for each pair is 88 m. 
The Daya Bay cores are separated from the Ling Ao 
cores by about 1100 m, while the Ling Ao-II cores are 
around 500 m away from the Ling Ao cores. 

Table 1. Vertical overburden, muon rate 
and average muon energy < _B M > of the three 
EHs. 





Overburden (m.w.e) 


i? M (Hz/m 2 ) 


< > (GeV) 


EH1 


250 


1.27 


57 


EH2 


265 


0.95 


58 


EH3 


860 


0.056 


137 



Three underground experimental halls (EHs) are 
connected with horizontal tunnels. For this analy- 
sis, two antineutrino detectors (ADs) were located 
in EH1, one in EH2, and three near the oscillation 
maximum in EH3 (the far hall). The overburden in 
equivalent meters of water (m.w.e.), simulated muon 
rate and average muon energy are listed in Table [T] 
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The distances from the six ADs to the six cores 
are listed in Table [2] All distances have been sur- 
veyed with the Global Positioning System (GPS) and 
with modern theodolites utilizing two major control 
networks built over several months. The network sur- 
veyed using GPS is within the campus of the power 
plant but outside of the tunnel. The other network 
is inside the tunnel system, surveyed using Total Sta- 
tion, an electronic/optical instrument widely used in 
modern surveying. The double traverse survey net- 
work was laid down in a closed ring in the 7-m wide 
tunnels. The Total Station survey included the power 
plant campus to link the two control networks. The 
survey from the anchors at the entrance of each ex- 
perimental hall to each AD was completed during the 
installation of each AD using a laser tracker. The co- 
ordinates of the AD center were further deduced us- 
ing the AD survey data collected during AD assembly. 
The coordinates of the geometrical center of the reac- 
tor cores were provided by the power plant relative to 
four anchor points outside of each nuclear island. The 
survey data were processed independently by three 
groups with different software. The uncertainty of 
the baselines was determined to be 28 mm as reported 
in Ref. Recently another closed traverse survey 
was completed utilizing a different tunnel entrance 
and the top of the mountain. The largest baseline 
difference between the two surveys is 4 mm and the 
uncertainty in the baselines has been reduced to 18 
mm. The uncertainty has seven significant contribu- 
tions, the largest being 12.6 mm due to the precision 
of the GPS survey. The second largest is 9.1 mm 
due to fitting uncertainties associated with the link- 
ing of the GPS and the Total Station networks. When 
combined with the uncertainties of the fission gravity 
center (described in Sec. 6), the baseline uncertainties 
were found to make a negligible contribution to the 
oscillation uncertainties. 

Table 2. Baselines from antineutrino detectors 
AD1-6 to reactors Dl, D2, and Ll-4 in meters. 





Dl 


D2 


LI 


L2 


L3 


L4 


AD1 


362 


372 


903 


817 


1354 


1265 


AD2 


358 


368 


903 


817 


1354 


1266 


AD3 


1332 


1358 


468 


490 


558 


499 


AD4 


1920 


1894 


1533 


1534 


1551 


1525 


AD5 


1918 


1892 


1535 


1535 


1555 


1528 


AD6 


1925 


1900 


1539 


1539 


1556 


1530 



uid scintillator (Gd-LS) [18j, [19|. The coincidence of 
the prompt scintillation from the e + and the delayed 
neutron capture on Gd provides a distinctive F e sig- 
nature. The positron carries almost all of the kinetic 
energy of the antineutrino, thus the positron energy 
deposited in the liquid scintillator is highly correlated 
with the antineutrino energy. The neutron thermal- 
izes before being captured on either a proton or a 
gadolinium nucleus with a mean capture time of ^30 
/is in Gd-LS with 0.1% Gd by weight. When a neu- 
tron is captured on Gd, it releases several gamma-rays 
with a total energy of ~8 MeV, and is thus easily dis- 
tinguished from the background coming from natural 
radioactivity. Only neutrons that captured on Gd 
were selected as the delayed signal of a antineutrino 
event in this analysis. 

Each AD has three nested cylindrical volumes sep- 
arated by concentric acrylic vessels [i^] as shown in 
Fig. [3J The innermost volume holds 20 t of Gd-LS 
with 0.1% Gd by weight and serves as the antineu- 
trino target. The middle volume is called the gamma 
catcher and is filled with 20 t of un-doped liquid scin- 
tillator (LS) for detecting gamma-rays that escape 
the target volume. The outer volume contains 37 t 
of mineral oil (MO) to provide optical homogeneity 
and to shield the inner volumes from radiation origi- 
nating, for example, from the photo-multiplier tubes 
(PMTs) or the stainless steel vessel (SSV) . There are 
192 20-cm PMTs (Hamamatsu R5912) installed along 
the circumference of the SSV and within the mineral 
oil volume, in 24 columns and 8 rings. To improve 
optical uniformity, the PMTs are recessed in a 3-mm 
thick black acrylic cylindrical shield located at the 
equator of the PMT bulb. 




2.2 Antineutrino Detectors 



The is e s are detected via the inverse /3-decay (IBD) 



reaction, v e + p 



■ n, in gadolinium-doped liq- 



Fig. 2. Schematic diagram of the Daya Bay de- 
tectors. 

Three automated calibration units (ACU-A, 
ACU-B, and ACU-C) are mounted on the top of each 
SSV as shown in Fig.O Each ACU is equipped with a 
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LED, a 68 Ge source, and a combined source of 241 Am- 
13 C and 60 Co. The Am-C source generates neutrons 
at a rate of -0.5 Hz. The rates of the 60 Co and 68 Ge 
sources are about 100 Hz and 15 Hz, respectively. 
Since the AD is fully submerged in water, the ACUs 
are operated remotely. The sources can be deployed 
to better than 0.5 cm along a vertical line down to the 
bottom of the acrylic vessels. When not in use, the 
LED and sources are retracted into the ACUs that 
also serve as shielding for the sources. 

2.3 Muon System 

The muon detection system consists of a resistive 
plate chamber (RPC) tracker and a high-purity ac- 
tive water shield. The water shield consists of two 
optically separated regions known as the inner (IWS) 
and outer (OWS) water shields. There are 121 (160) 
PMTs installed in the IWS and 167 (224) PMTs in 
the OWS in each near (far) hall. Each region op- 
erates as an independent water Cherenkov detector. 
The muon detection efficiency is 99.7% and 97% for 
the IWS and OWS, respectively [lj|. In addition to 
detecting muons that can produce spallation neutrons 
or other cosmogenic backgrounds in the ADs, the pool 
moderates neutrons and attenuates gamma rays pro- 
duced in the rock or other structural materials in and 
around the experimental hall. At least 2.5 m of wa- 
ter surrounds the ADs in every direction. Each pool 
is outfitted with a light-tight cover overlaying a dry- 
nitrogen atmosphere. 

Each water pool is covered with an array of RPC 
modules 2l|, [22j . The 2 m x 2 m modules are lay- 
ered on a steel frame to minimize dead areas. The 
assembly is mounted on rails and can be retracted to 
provide access to the water pool. There are four lay- 
ers of bare RPCs inside each module, with one layer 
of readout strips associated with each layer of bare 
RPCs. The strips have a "switchback" design with 
an effective width of 25 cm, and are stacked in alter- 
nating orientations providing a spatial resolution of 
~8 cm. 

2.4 Trigger and Readout 

Each detector unit (AD, IWS, OWS, and RPC) is 
read out with a separate VME crate. All PMT read- 
out crates are physically identical, differing only in 
the number of instrumented readout channels. The 
front-end electronics board (FEE) receives raw sig- 
nals from up to sixteen PMTs, sums the charge from 
all input channels, identifies over-threshold channels, 
records their timing information, and measures the 
charge of each over-threshold pulse with a 40 MHz 



sampling rate 23J. The FEE in turn sends the num- 
ber of channels over threshold and the integrated 
charge to the trigger system. When a trigger is issued, 
the FEE reads out the charge and timing information 
within 1 /us for each over-threshold channel, as well 
as the average ADC value over a 100 ns time- window 
immediately preceding the over-threshold condition 
(preADC). 

Triggers are primarily created internally within 
each PMT readout crate based on the number of over- 
threshold channels (NHIT) as well as the summed 
charge (E-Sum) from each FEE [24]. The system is 
also capable of accepting external trigger requests, for 
example, from the calibration system. The trigger 
system blocks triggers when either the trigger data- 
buffer or a FEE data-buffer is nearly full. The num- 
ber of blocked triggers is recorded and read out for 
calculating the dead time offline. 

3 Data Characteristics, Calibration 
and Modelling 

3.1 Data set 

The data used in this analysis were collected 
from December 24, 2011 through May 11, 2012. Ta- 
ble [3] summarizes the experimental livetime for each 
hall. Total data acquisition (DAQ) time measures the 
number of hours that the DAQ was collecting data, 
with about 2% of the DAQ time devoted to detector 
calibration. Standard data running (Physics Data or 
Physics DAQ time) accounted for more than 93% of 
the calendar time. We further rejected about 60 hours 
of physics data from each hall due to excessive coher- 
ent electromagnetic pickup, PMT high voltage (HV) 
trips, electronic or DAQ problems, or requirements of 
simultaneous operation in all three halls. The result- 
ing data set (Good run data or Good run time) were 
used for analysis. 



Table 3. 
hours. 



Summary of experimental livetime in 





EHl 


EH2 


EH3 


Total calendar time 


3322.1 


3322.1 


3322.1 


Total DAQ time 


3195.4 


3179.5 


3171.6 


Physics DAQ time 


3117.9 


3122.0 


3093.6 


Good run time 


3061.1 


3057.1 


3030.5 



The detector halls operated independently with 
a common centralized clock and GPS timing sys- 
tem. The analysis presented here required simulta- 
neous operation of all three detector halls, to mini- 
mize systematic effects associated with potential re- 
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actor power excursions. Simultaneous operation was 
defined as Physics Data within a given hour existing 
for all three detector halls. The data samples used in 
this analysis differed by 1% in time for the three halls. 
A more rigorous requirement that demands synchro- 
nization among the three halls on the scale of seconds 
was tested with no change to the reported results. 

3.2 Triggered Detector Rates 

Triggers were formed based either on the number 
of PMTs with signals above a ^0.25 photoelectron 
(p.e.) threshold (NHIT triggers), or the charge sum of 
the PMTs (E-Sum triggers). AD triggers with NHIT 
> 45 or E-Sum > 65 p.e. correspond to an event en- 
ergy threshold of ~0.4 MeV [3]. The corresponding 
trigger rate per AD was < 280 Hz with a negligible 
trigger inefficiency for IBD candidates. 

The I7 e candidates were selected in the offline anal- 
ysis using the coincidence of a prompt signal from the 
e + and a delayed signal due to neutron capture on Gd. 
A prompt-type (delayed-type) signal was defined as 
an event with energy in the range of 0.7-12 MeV (6- 
12 MeV). The rates of prompt- type and delayed-type 
singles that are separated in time by at least 200 as 
from any additional signals with an energy > 0.7 MeV 
were of particular interest for background studies and 
detector stability monitoring. They are shown in 
Fig. |31 A veto was applied to reject events within -2 
to 200 ^s relative to a muon (defined in Sec. 4.1). The 
data were corrected for the corresponding inefficien- 
cies. These rates were used to estimate the accidental 
background rate as described in Sec. 5.1. 



♦ ADl -©- AD4 
*AD2 -&AD5 
-*-AD3 *AD6 




Jan 21 



Feb 20 



Mar 21 



Apr 20 



Date 



Fig. 3. Singles rates for the six ADs. The top 
panel shows the prompt candidates and the 
bottom panel shows the delayed candidates. 

The observed rate of low energy signals decreased 
with time. The detectors in EH1 initiated data tak- 



ing on Aug. 15, 2011 and the AD in EH2 started on 
Nov. 5, 2011. As such, these detectors (AD1-3) had 
reached a steady state by December 24, 2011, while 
the rates in AD4-6 in EH3 exhibited decaying behav- 
ior, as shown in Fig. [3] 

The muon rates in the water Cherenkov detectors 
(IWS and OWS) were closely monitored, as shown 
in Fig. 0J IWS and OWS events were selected with 
NHIT > 12. The event rates were different for the 
three halls due to differing muon rates in each hall 
and different sizes of the far hall and the near halls. 



220 
200 
180 
160 



140 
120 
100 
80 
60 
40 
20 



i.iii.i.Tii.ii j-ii lifc.'iii' ii I'li-r- lini ii 



-iniiiium 



-•-EH1 IWS 
-&-EHI OWS 



- EH2 IWS -•- EH 3 IWS 

- EH2 OWS -b- EH 3 OWS 



J I I I I I , I I I , I I I I I , ■_ 



Apr 20 



Date 

Fig. 4. Muon rates in the inner (IWS) and 
outer water shield (OWS) in the three exper- 
imental halls. 

3.3 Instrumental Backgrounds 

A small number of AD PMTs spontaneously emit 
light, due to discharge within the base. These instru- 
mental backgrounds are referred to as flasher events. 
For Daya Bay, the reconstructed energy of such events 
covers a wide range, from sub-MeV to 100 MeV. Two 
features were typically observed when a PMT flashed: 
the observed charge fraction for a given PMT was 
very high, and PMTs on the opposite side of the AD 
saw large fraction of light from the flashed PMT. The 
charge pattern of a typical flasher event is shown in 
Fig. El 

To reject flasher events, two variables, named 
MaxQ and Quad, were created based on the distinc- 
tive charge pattern. MaxQ is the largest fraction of 
the total detected charge seen by a single PMT (the 
"hottest" PMT). There are twenty-four columns of 
PMTs in an AD that can be divided into four quad- 
rants. With the hottest PMT centered in the first 
quadrant, Quad was defined as Q3/ (Q2 + Q4), where 
Qi is the charge sum of the PMTs in the i-th quad- 
rant. A flasher event identification variable (FID) was 
constructed based on MaxQ and Quad: 



FID = log 10 [(MaccQ/0.45) 2 + (Quad) 2 



(2) 
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Fig. [6] shows the discrimination of flasher events for 
the delayed signal of the IBD candidates. The distri- 
butions for all six ADs agree well for IBD candidates 
(FID<0); however, there is some variation for flasher 
candidates (FID>0). For the IBD analysis as well 
as other analyses, the rejection of flasher events was 
done at the beginning of the data reduction. 




Fig. 5. Topology of a typical flasher event. 
Such events are distinctive, characterized by a 
single channel with substantially more charge 
than in surrounding PMTs, as well as exces- 
sive charge on the opposite side of the AD. 




-1.5 



-0.5 



0.5 1 1.5 2 
FID of delayed candidates 

Fig. 6. Discrimination of flasher events (FID 
>0) and IBD delayed signals (FID <0). The 
delayed signals of IBDs have the same distri- 
bution for all six ADs while the flashers are 
different. The FID < distributions have been 
scaled to equal area. 

The discrimination power decreases for low en- 
ergy events or events very close to PMTs. For the 
rejected events with FID ~ 0, we studied the charge 
pattern, the energy distribution, the capture time, 
and the distance between the vertices of the prompt 
and the delayed signals, and found that some were 
consistent with real IBD events. By counting such 
events, the inefficiency of the IBD selection due to the 
flasher rejection was estimated to be 0.02% with an 
uncorrelated uncertainty of 0.01%. The background 



contamination of selected IBD candidates was evalu- 
ated to be < 10 -4 . Furthermore, such contamination 
was counted as accidental background (Sec. 5.1) and 
was subtracted. Special runs were conducted with re- 
duced high voltage for selected PMTs to cross-check 
the identified PMTs that exhibited flashing. Due to 
the high efficiency of the FID, all AD PMTs were 
kept in operation, including those identified as flash- 
ing PMTs. 

3.4 Energy Reconstruction 

In general, the energy response of the AD can de- 
pend on time, position in the fiducial volume (non- 
uniformity), particle species, and their energies (non- 
linearity). The goal of energy reconstruction was to 
correct for these dependencies in order to minimize 
the uncertainties in the AD energy scale. To achieve 
this goal, each AD was calibrated using LEDs, 68 Ge, 
and 241 Am- 13 C/ 60 Co sources. LEDs were utilized for 
PMT gain calibration, while the energy calibration 
parameter (p.e. per MeV) was determined with a 60 Co 
source deployed at the detector center. The sources 
were deployed once per week to check and correct 
for any time dependence. Occasionally a PMT's out- 
put was noisy and was temporarily turned off dur- 
ing physics data taking. The energy calibration cor- 
rected such situations. The energy calibration param- 
eter for each AD is shown in Fig. [7] as a function of 
time. The small jumps correspond to the temporary 
turn-off of noisy PMTs. The energy resolution was 
(7.5/ x / J B(MeV) + 0.9)% for all 6 ADs. 
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Fig. 7. Calibration parameter versus time for 
each AD. 

A scan along the vertical axis utilizing the 60 Co 
source from each of the three ACUs was used to ob- 
tain non-uniformity correction functions. The non- 
uniformity was also studied with spallation neutrons 
generated by cosmic muons, and alphas produced by 
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natural radioactivity present in the liquid scintillator. 
The neutron energy scale was set by comparing 60 Co 
events with neutron capture on Gd events from the 
Am-C source at the detector center. Additional de- 
tails of energy calibration, reconstruction, and vertex 
reconstruction can be found in Ref. 



3.5 Detector Simulation 



13- 



The AD energy scale uncertainty was studied by 
comparing the energy peaks in all ADs using neutron 
capture on gadolinium from IBD and muon spallation 
products, alphas from Polonium decay in the Gd-LS, 
and each of the calibration sources. Asymmetries of 
the six ADs' response are shown in Fig. El For each 
type of event, we defined the Asymmetry as 



Asymmetryi — 



(3) 



where E t is the fitted mean energy of the studied type 
of event of the i-th AD. The energy scale uncertainty 
was set at 0.5% in Ref. [iff based on extensive side- 
by-side studies of AD1 and AD2. Extending this to 
six ADs, Asymmetries for all types of events in all the 
ADs fall within a band of 0.5%. As such, we kept the 
same uncertainty, uncorrelated among ADs. 



>, 0.01 
| 0.008 
| 0.006 
< 0.004 
0.002 


-0.002 
-0.004 



,5 Ge 
'°Co 



Po a 
4 Poa 



a Am-Cn-Gd * 2U Poa 



ir IBD n-Gd 
★ Spa n-Gd 



ft" 



5 6 
AD Number 



b 0.006 p 
| 0.004 §- 
1 0.002 i- 

I* 

^ -0.002 r 
-0.004 r- 
-0.006 5- 







6 8 10 

Reconstructed Energy (MeV) 



Fig. 8. Asymmetries in energy response for all 
six ADs. The sources 68 Ge, 60 Co, and Am- 
C were deployed at the detector center. The 
60 Co data was used for energy calibration. 
The alpha particles from Polonium decay and 
neutron capture on gadolinium of IBD and 
spallation neutrons were uniformly distributed 
within each detector. Differences between 
these sources are due to spatial non-uniformity 
of the detector response. The same set of data 
points are shown in the lower panel as a func- 
tion of energy, which demonstrates that all six 
ADs have similar energy non-linearity. 



A Geant4 [25j based computer simulation (Monte 
Carlo, MC) of the detectors and readout electron- 
ics was used to study the detector response. It con- 
sisted of five components: kinematic generator, detec- 
tor simulation, electronics simulation, trigger simula- 
tion and readout simulation. The MC was carefully 
tuned to match observed detector distributions, such 
as PMT timing, charge response, and energy non- 
linearity. 

The antineutrino generator read from a database 
that stored the reactor antineutrino spectra from each 
core at each detector. The database was binned in 
daily increments and accounted for fuel evolution. 
The flux was scaled later based on the actual reactor 
power. The cosmic muons in the underground labo- 
ratory were simulated using a digitized to pog raphic 
map of the site and Muon Simulation Code [26J (MU- 
SIC), which calculated the energy loss and multiple 
scattering due to the rock overburden. The muon 
generator for Geant4 read randomly from a library 
of muon events generated with MUSIC. The software 
generators for the calibration sources and the simula- 
tion of the decay sequences for natural radionuclides 
found in our detectors were customized based on data 
from the ENDF database 



27] 



All physical processes in Geant4 relevant to the 
Daya Bay simulation were validated. In the valida- 
tion process, we found that the gamma spectra of 
neutron capture and muon capture on many nuclei 
were incorrectly modeled. Since a systematic correc- 
tion was complex, we implemented corrections on a 
case by case basis. The most important one was the 
neutron capture on gadolinium where we used a cus- 
tomized module based on the measured gamma spec- 
trum [jjj ]. Furthermore, the simulation of thermal 
neutron scattering was improved by considering the 
molecular binding energy of the scattering nuclei. 

The gadolinium and other elemental concentra- 
tions of the liquid scintillator were measured and in- 
corporated into the MC. All relevant optical prop- 
erties of the detector components were derived from 
measurements, including the attenuation lengths and 
refractive indices of all liquids as well as the acrylic 
components, time constants and photon emission 
spectra of Gd-LS, LS, and mineral oil, and the reflec- 
tivity of the reflectors as well as other detector ma- 
terials. Photon absorption and re-emission processes 
in liquid scintillator were modeled based on measure- 
ments in order to properly simulate the propagation 
of optical photons and contributions from Cherenkov 
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process. 

The details of the electronics simulation can be 
found in Ref. [29(. Using the timing and number 
of p.e. generated in PMTs, an analog signal pulse 
for each PMT was generated and tracked through 
the digitization process, taking into account the non- 
linearity, dark rate, pre-pulsing, after-pulsing and 
ringing of the waveform. The simulated analog pulse 
was then used as input to a trigger system simulation 
for each sub-detector. 

4 Event Selection 

4.1 IBD Selection 

Two conditions were implemented prior to the 
IBD selection. First, flasher events were rejected 
(Sec. 3.3). Second, all AD triggers within a (-2 /is, 
200 (is) time-window with respect to a water shield 
muon candidate ((iws) were rejected, where a (i W s 
was defined as any signal with NHIT > 12 in ei- 
ther the inner or outer water shield. This allowed for 
the removal of most of the superfluous triggers that 
followed a muon, as well as triggers associated with 
muon-induced spallation products. The veto time- 
window was extended to 2 /is earlier than the muon 
to avoid time alignment issue among different detec- 
tors. Events in an AD within ±2 /is of a (iws with 
energy >20 MeV or >2.5 GeV were classified as AD 
muons ((Iad) or showering muons (/i S h), respectively. 
Longer veto windows were applied for such events to 
further reject cosmogenic backgrounds. 

The energy of the prompt and delayed candidates 
were required to satisfy 0.7 MeV < E v < 12.0 MeV 
and 6.0 MeV < E d < 12.0 MeV, respectively, and 
At = t d — t p must have satisfied a 1 < At < 200 /is 
coincidence, where t p and t d are the times of the 
prompt and delayed signals. A multiplicity cut re- 
quired no additional candidate with E > 0.7 MeV in 
the interval 200 /is before t p , 200 (jls after t d , or be- 
tween t p and t d . The prompt-delayed pair was vetoed 
if the delayed candidate satisfied any of the condi- 
tions — 2 (j,s < t d — t flws <600 /is (water shield muon), 
< t d -t MAD < 1000 (is (AD muon), or < t d —t^ h <l s 
(AD showering muon). The prompt energy, delayed 
energy and capture-time distributions for data and 
MC are shown in Figs. [9] - [TT1 respectively. 

The data are generally in good agreement with the 
MC. The apparent difference between data and MC 
in the prompt energy spectrum in Fig. [9] is primarily 
due to nonlinearity of the detector response. Since all 
ADs had similar nonlinearity (as shown in the bot- 



tom panel of Fig. [8J, and the energy selection cuts 
cover a larger range than the actual distribution, the 
discrepancies between MC and data introduced neg- 
ligible uncertainties to the rate analysis. Therefore, 
this nonlinearity correction was not implemented in 
this analysis. 
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Fig. 9. Prompt energy spectrum from AD1. 
IBD selection required 0.7 < E p < 12.0 MeV. 
The spectrum of accidental backgrounds, de- 
termined from the distribution of all prompt- 
type signals, was subtracted. 
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Fig. 10. Delayed energy spectrum from AD1. 
IBD selection required 6.0 < E d < 12.0 MeV. 
The spectrum of accidental backgrounds, de- 
termined from the distribution of all delayed- 
type signals, was subtracted. 

4.2 Efficiencies and Uncertainties 

For a relative measurement, the absolute efficien- 
cies and correlated uncertainties do not factor into 
the error budget. In that regard, only the relative ef- 
ficiencies and uncorrelated uncertainties matter. Ex- 
traction of absolute efficiencies and correlated uncer- 
tainties was done in part to better understand our 
detector, and was a natural consequence of evaluat- 
ing the uncorrelated uncertainties. Absolute efficien- 
cies associated with the prompt energy, delayed en- 
ergy, capture time, Gd-capture fraction, and spill-in 
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effects were evaluated with the Monte Carlo. Efficien- 
cies associated with the muon veto, multiplicity cut, 
and livetime were evaluated using data. In general, 
the uncorrelated uncertainties were not dependent on 
the details of our simulation. 
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Fig. 11. Neutron capture time from AD1. IBD 
selection required 1 < td — t p < 200 [is. In or- 
der to compare data with MC, a cut on the 
prompt energy (E p > 3 MeV) was applied to 
suppress accidental backgrounds. A zoomed- 
in plot for 1 < td — t p < 30 /is is shown in the 
inset. 

Table [4] summarizes the absolute efficiencies and 
the systematic uncertainties. The uncertainties of the 
absolute efficiencies were correlated among the ADs. 
No relative efficiency, except the muon veto efficiency 
and the average multiplicity cut efficiency Z m , were 
corrected. All differences between the functionally 
identical ADs were taken as uncorrelated uncertain- 
ties. 

Table 4. Summary of absolute efficiencies, and 
correlated and uncorrelated systematic uncer- 
tainties. For our relative measurement, the 
absolute efficiencies as well as the correlated 
uncertainties effectively cancel. Only the un- 
correlated uncertainties contribute to the final 
error in our relative measurement. 





Efficiency 


Correlated 


Uncorrelated 


Target protons 




0.47% 


0.03% 


Flasher cut 


99.98% 


0.01% 


0.01% 


Delayed energy cut 


90.9% 


0.6% 


0.12% 


Prompt energy cut 


99.88% 


0.10% 


0.01% 


Multiplicity cut 




0.02% 


<0.01% 


Capture time cut 


98.6% 


0.12% 


0.01% 


Gd capture fraction 


83.8% 


0.8% 


<0.1% 


Spill-in 


105.0% 


1.5% 


0.02% 


Livetime 


100.0% 


0.002% 


<0.01% 



Combined 78.8% 1.9% 0.2% 



The absolute efficiency of the prompt energy cut 
(0.7 < E p < 12.0 MeV) was determined to be 99.88%. 
The energy spectrum is shown in Fig. |9l Inefficiency 
was mainly caused by interactions inside the inner 
acrylic vessel, indicated by the vertex distribution 
of the rejected prompt signal below 0.7 MeV. While 
the uncertainty in the energy scale is below 0.5% for 
events at the detector center or uniformly distributed 
in the target volume, it is larger for events at the 
edge. The asymmetries of the energy among ADs 
could be as large as 2% for events at the radius of 
ACU-C, studied with 60 Co source 0. Taking 2% 
uncertainty in the energy scale for events near the in- 
ner acrylic vessel, the uncorrelated uncertainty of the 
efficiency due to the prompt energy cut was evalu- 
ated to be 0.01%. Given that the positron threshold 
was calibrated with the 68 Ge source, the uncertainty 
of this absolute efficiency comes from the difference 
of non-linearity and non-uniformity between the data 
and MC. The correlated uncertainty was estimated 
to be 0.1%. 

The absolute efficiency of the delayed energy cut 
(6.0 < E d < 12.0 MeV) was determined to be 90.9%. 
As shown in Fig. [TU1 the fraction of events in the 6-7 
MeV region was 5.3% of that in 6-12 MeV for MC. 
For selected IBD data, this fraction was 5.6%. As- 
suming the same relative difference between MC and 
data in the 0-6 MeV region, the difference of absolute 
efficiency between the MC and data was evaluated to 
be 0.6%, which is taken as the correlated uncertainty. 
By varying the cut at 6 MeV and counting the num- 
ber of events in the selected sample, we found that the 
0.5% asymmetry of the energy scale in ADs, shown 
in Fig. [SI leads to a 0.12% uncorrelated efficiency un- 
certainty. The low energy tail around 6 MeV is rel- 
atively flat and the MC and data agree well. Both 
MC and data studies yield the same uncorrelated ef- 
ficiency uncertainty. 

The spill-in enhancement resulted when neutrons 
from IBD interactions outside the target volume were 
captured by a Gd nucleus in the target volume. It was 
defined as the ratio of all IBD interactions that lead to 
a neutron capture on Gd to IBD interactions within 
the target volume leading to a neutron capture on Gd. 
From MC, it was evaluated to be 105.0%. By model- 
ing the relative difference in acrylic vessel thickness, 
acrylic density and liquid density in MC, the relative 
uncertainty of the spill-in efficiency was evaluated to 
be 0.02%. The correlated uncertainty of the spill-in 
efficiency was evaluated with MC. The modeling of 
molecular binding energy of the scattering nuclei has 
a large impact on the simulation of thermal neutron 
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scattering, and thus on the absolute spill-in efficiency. 
The thermal neutron scattering process is correlated 
with the neutron capture time. The agreement be- 
tween data and MC is shown in the inset of Fig. [TT] 
By comparing the results of simulation with two dif- 
ferent models of molecular binding energy as well as 
without binding energy, we conservatively estimated 
the correlated uncertainty of the spill-in efficiency to 
be 1.5%. 

The Gd capture fraction was defined as the ra- 
tio of the number of Gd capture events produced 
by IBD reactions to all IBD reactions in the Gd-LS. 
It was evaluated to be 83.8%. The spill-out deficit, 
^2.2% by comparing the Gd capture fraction of the 
Am-C neutron source at the detector center and IBD 
events in MC, was included in the absolute Gd cap- 
ture fraction. Spill-out is analogous to spill-in, except 
that IBD neutrons produced within the target volume 
were captured outside the target volume. By measur- 
ing the difference in the neutron capture time of each 
AD, the relative Gd-concentration variation was con- 
strained and the Gd capture fraction variation was 
determined to be within 0.1%. By comparing Am-C 
source data with MC, as well as spallation neutrons, 
the correlated uncertainty on Gd capture fraction was 
estimated to be 0.8%. 

The efficiency of the capture time cut (1 < At < 
200 /is) was evaluated to be 98.6% with 0.2% of events 
with At < 1 (is and 1.2% events with At > 200 /is. 
The correlated efficiency uncertainty was evaluated 
to be 0.12%, according to the difference in the mea- 
sured capture time between Am-C data and MC. 
The uncorrelated uncertainty comes from the Gd- 
concentration variation and possible trigger time- 
walk effect, and it was evaluated to be 0.01%. 

The muon veto efficiencies were determined using 
data. For each type of muon candidate (/i ff s, [J* ad 
and /J-sh), the start and end time of the veto window 
were well defined. Overlapping veto windows were 
merged to avoid double counting. As a result each 
livetime window was precisely calculated as the un- 
vetoed time interval between two isolated veto win- 
dows. The total livetime was obtained by summing 
all the individual livetime windows. The muon veto 
efficiency e M was defined as the fraction of the live- 
time after a muon veto in the total DAQ livetime. 
For each experimental hall the muon rates were sta- 
ble as shown in Fig. The muon veto efficiencies 
differed due to different muon candidate rates. 

The multiplicity cut required no additional 
> 0.7 MeV signals (singles) in the time range from 
200 /is before the prompt signal to 200 /is after the 



delayed signal. The singles rate R a can be taken 
as the rate of prompt-type signals shown in Fig. [3l 
The efficiency of the multiplicity cut is a product 
of three components. The probability of no singles 
in the 200 /is before the prompt signal is given by 
cxp(— i?2oo), where i? 2 oo = Rs • 200/is. The prob- 
ability of no singles between the prompt and de- 
layed signal is given as / ^ exp(—R 3 t)f(t)dt, where 
f(t) is the probability density function of the cap- 
ture time of neutron on Gd, and can be simplified as 
1 — R s t C np + O(10 -5 ), where i cap is the mean neutron 
capture time in 200 /is. The average of the mean cap- 
ture time of the six ADs was 33.46 /is, obtained from 
data. The uncorrelated uncertainty was determined 
by the difference of the mean capture times among 
ADs. The probability of no singles in 200 /is after 
the delayed signal must be calculated for two cases 
since the window may be truncated by an AD muon 
that would obscure any potential single. If the single 
livetime window was T 3 < 200/is, the efficiency was 



1 



R s T a 

and if T s > 200/is, the efficiency was 

1 



(1 t^) e - R 200 



(i- 



-R-200 



Because the second term depends on the length of the 
single livetime window, the multiplicity cut efficiency 
must be calculated for every single livetime window. 
As a consequence, the muon veto efficiency and the 
multiplicity were coupled. The combined efficiency is 



DAQ ! 



(4) 



where is the multiplicity cut efficiency in the i-th 
single livetime l 7 *, and Tb A Q is the analyzed good run 
time. The muon veto efficiency e M and the average 
multiplicity cut efficiency e m calculated with Eq. 2] 
arc listed in Table 5 and corrected for each AD. 

The target mass uncertainty was discussed exten- 
sively in Ref. [laj . The correlated uncertainty 0.47% 
largely comes from the hydrogen-carbon ratio of the 
target liquid, which is canceled out in the near-far 
relative measurement by using the same batch of Gd- 
LS. The time variation of the target mass, e.g. due 
to temperature variation, is monitored by the liquid 
level with several independent sensors in the overflow 
tanks on the top of the AD lid [30| . The variation 
of the target mass for the analyzed data set is shown 
in Fig. [12j The ±0.02% range is the target mass un- 
certainty evaluated during filling [15j . To accurately 
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evaluate the mass of Gd-LS transferred into detec- 
tors, a 20-t filling tank was equipped with load cells 
to measure the mass of the filling tank before and af- 
ter filling. The above uncertainty is dominated by the 
load cell drift during the filling operation. As such, 
the uncorrelated uncertainty is set to be 0.03%. 
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Fig. 12. Target mass variation for each AD over 
the analyzed time period. The vertical double 
arrow indicates the total uncorrelated uncer- 
tainty in the target mass evaluated during fill- 
ing. 

5 Backgrounds 

5.1 Accidental Backgrounds 

The accidental background is defined as any pair 
of otherwise uncorrelated signals that happen to sat- 
isfy the IBD selection criteria. For any given sig- 
nal with an observed energy between 6 and 12 McV 
(delayed- type signal), the probability of forming an 
accidental background is the product of two compo- 
nents, the probability of a prompt-type signal within 
1-200 fis before the delayed-type signal, 1— exp(— R- 
199/is), and the probability of no singles within 200/is 
before the prompt-type signal and 200 /is after the 
delayed-type signal, exp(— R-400fj,s). R is the rate of 
prompt-type singles. Since the rate of prompt-type 
and delayed-type singles changed over time, the ac- 
cidental background was calculated every four hours 
and summed as follows: 



iVacc. bkg .=5>,e-*> 



■400/js 



(I" 



(5) 



where N t and Ri are the number of delayed-type and 
prompt-type singles rates in the i-th four-hour period, 
respectively. The statistical uncertainty was domi- 
nated by Ni, and was approximated as 

A ' acc. bkg. 



r Ar(stat.) 
0Jv acc.bkg. 



(6) 



The expected rates of accidental backgrounds are 
listed in Table [H after correcting for the muon veto ef- 
ficiency and the multiplicity cut efficiency in the IBD 
selection. 

An alternate method to determine the acciden- 
tal backgrounds, the off-window method, was de- 
veloped. By definition, the accidental background 
within the IBD coincidence time window (1 fis < At < 
200 fis) should be the same as in any other window 
(t oH + 1 fis < At < t ott + 200 fis), where t oS is an ar- 
bitrary time offset. If t° a is large enough to avoid 
real correlated events (such as for IBD, fast neutron 
(Sec. 5.2), and 9 Li/ 8 He decay (Sec. 5.3)), the acci- 
dental backgrounds can be estimated by counting the 
coincidences in the off-window. To reduce the statisti- 
cal uncertainty, multiple non-overlapping off- windows 
were examined. The mean number of selected coinci- 
dences in these off- windows was taken as the expected 
accidental background. The relative differences be- 
tween the results from the off-window method and 
the calculations using Eq. [5] were consistent given the 
statistical uncertainties for all six ADs. 

The accidental background was also validated by 
comparing the distributions of distance between the 
reconstructed vertices for the prompt and delayed sig- 
nals of the IBD candidates and accidentals selected 
by the off-window method, as shown in Fig. 1131 The 
prompt and delayed vertices of accidentals were un- 
correlated, thus giving a broad distribution, while the 
two vertices for IBD events were correlated, giving a 
distribution peaked at short distance. For distances 
greater than 2 m, the IBD candidate and off-window 
distributions agree well. 

5.2 Fast Neutron Backgrounds 

Energetic neutrons created by cosmic rays enter- 
ing an AD could mimic IBD by recoiling off a proton 
before being captured on Gd. Since the visible en- 
ergy of the recoil proton ranged well past that of the 
IBD events (up to several hundred MeV as shown 
in Fig. I14p . we estimated the number of fast- neutron 
background events in the IBD sample by extrapolat- 
ing the prompt energy (E p ) distribution between 12 
and 100 MeV down to 0.7 MeV. Two different ex- 
trapolation methods were used. By assuming the 
recoil proton energy spectrum follows a flat distri- 
bution, the mean number of events per energy bin 
of the distribution from 12 to 100 MeV was used to 
estimate the number of fast-neutron events between 
0.7 and 12 McV. Alternately, the data from 12 to 
100 MeV were fit with a first-order polynomial func- 
tion (f(E) = a + bE). The best-fit parameters were 
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Table 5. Summary of signal and background. The background and IBD rates have been corrected for the 
muon veto efficiency e M and the average multiplicity cut efficiency e m . 





AD1 


AD2 


AD3 


AD4 


AD5 


AD6 


IBD candidates 


69121 


69714 


66473 


9788 


9669 


9452 


Expected IBDs 


68613 


69595 


66402 


9922.9 


9940.2 


9837.7 


DAQ livetime (days) 


127.5470 


127.3763 




126.2646 






0.8231 


0.8198 


0.8576 


0.9813 


0.9813 


0.9810 




0.9738 


0.9742 


0.9753 


0.9737 


0.9734 


0.9732 


Accidentals (per day) 


9.73±0.10 


9.61±0.10 


7.55±0.08 


3.05 ±0.04 


3.04 ± 0.04 


2.93 ±0.03 


Fast-neutron (per day) 


0.77±0.24 


0.77±0.24 


0.58±0.33 


0.05±0.02 


0.05±0.02 


0.05±0.02 


9 Li/ 8 He (per AD per day) 


2.9±1.5 


2.0±1.1 




0.22±0.12 




Am-C correlated (per AD per day) 






0.2±0.2 






(a, n) background (per day) 


0.08±0.04 


0.07±0.04 


0.05±0.03 


0.04±0.02 


0.04±0.02 


0.04±0.02 


IBD rate (per day) 


662.47±3.00 


670.87±3.01 


613.53±2.69 


77.57±0.85 


76.62±0.85 


74.97±0.84 



used to estimate the number of fast-neutron events 
between 0.7 and 12 MeV. The fast neutron back- 
ground in the IBD sample was assigned to be equal 
to the mean value of the two extrapolation methods. 
The systematic error was determined from their dif- 
ferences and the fitting uncertainties. 



E 




1 2 3 4 5 6 

Prompt-delayed Distance(m) 

Fig. 13. Distance between the prompt sig- 
nal and delayed signal. The dots show the 
IBD candidates in data and the open cir- 
cles are accidental candidates selected with 
the off-window method, both in their abso- 
lute rates. The histogram shows the simulated 
IBD events, with rate normalized to data. 

As a check, we studied the fast neutrons associ- 
ated with tagged muons. The prompt energy of the 
fast neutron tagged by the IWS muon will be con- 
taminated if the muon clips the edge or corner of an 
AD. Furthermore, the fast neutron backgrounds in 
the IBD candidate pool mostly originated from OWS 
muons (defined as OWS PMT multiplicity > 12 and 
without an IWS trigger) or muons passing through 
rock, since the muon detection efficiency of the IWS 
was very high (99.7%). The fast neutrons tagged by 
the OWS muons or RPC-only muons (only detected 



by RPC) had a prompt energy spectrum similar to 
the fast neutron backgrounds in the IBD sample. Af- 
ter rejecting flasher events, we selected fast-neutron- 
like events by requiring exactly two signals within 
200 /is after an OWS muon or a RPC-only muon. 
The time interval and the energy selections of the 
prompt-delayed pair were the same as the IBD selec- 
tions, except that the prompt energy was relaxed to 
be 0.7 < E p < 100 MeV. We combined the sam- 
ples from EH1 and EH2 to create the fast neutron 
prompt energy spectrum shown in the inset of Fig ll4l 
The observed distribution validates our extrapolation 
method for estimating the fast neutron background. 
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Fig. 14. Prompt energy spectrum of IBD can- 
didates with the upper limit relaxed. The 
energy spectrum of the fast neutron back- 
grounds tagged by the OWS muons or RPC- 
only muons is shown in the inset. 

Two additional methods were used to provide fur- 
ther cross checks and estimates for the fast neutron 
background. These two methods are consistent with 
the result by extrapolating the IBD prompt energy 
spectrum within the assigned uncertainty. 
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First, a muon with a large PMT multiplicity in 
the IWS and OWS has an increased probability to 
produce a fast neutron in an AD, presumably since 
track length correlates with PMT multiplicity. Such 
a correlation has been observed in data. Also, muon 
detection inefficiency was associated with low PMT 
multiplicity. By extrapolating the fast neutron rate 
produced by muons with a sum of PMT multiplicities 
between 24 and 48 in the IWS and OWS to the range 
and 24, we were able to estimate the fast neutron 
background slipping into the IBD sample due to the 
inefficiency of the muon detection. 

Second, we collected different fast neutron sam- 
ples based on muons going through different detec- 
tor volumes (n^ ws : fast neutron from an IWS tagged 
muon; n° ws : fast neutron from an OWS muon; 
n™ ck : fast neutron from a muon going through nearby 
rock) and estimated these samples separately. The 
n r f ° ck was estimated by selecting RPC-only muons. 
MC simulation suggested that the fast neutron back- 
grounds tagged by RPC-only muons account for one- 
third of the rock neutron background. The fast neu- 
tron background (n/) is described as 

n, = n™ s (l - 6ws) + < ws (1 - £ows) + n r / ck (7) 

where £iws is the muon detection efficiency of the 
IWS and Cows is that of the OWS. 

5.3 9 Li/ 8 He Backgrounds 

The rate of correlated background from the /3-n 
cascade of the cosmogenic 9 Li/ 8 He decays was evalu- 
ated from the distribution of the time since the last 
muon, which can be described as [iH 

f® = T~ e~ t/Xa + t- ■ e ~ t/Xb + %^"* /T , (8) 

""V Aft 1 

where B a and B b are the number of /3-n events for 9 Li 
and 8 He, respectively. T is the mean time between 
muons, t- = ± + — and -f- = ^ + — with r a = 0.257 
s and r b — 0.172 s being the known decay time con- 
stants for 9 Li and 8 He, respectively. The muon rate 
i? M = 1 /T depends on the muon selection criteria. 

To reduce the impact of accidental backgrounds 
on our measurement of 9 Li and 8 He, we made the 
following modification to our IBD selection criteria: 

• 0.7 <E p < 12.0 MeV changed to 3.5 < E p < 12.0 
MeV. 

• 1 < At < 200/is changed to 1 < At < 100/zs. 

• fi sh veto time changed from 1 s to 1000 /is 



The measured 9 Li/ 8 He rate was corrected for the rel- 
ative efficiency with respect to the IBD selection cri- 
teria. Assuming that 9 Li was predominant over 8 He 
(as observed in a previous experiment Q and consis- 
tent with our observations), and based on the 9 Li (3 
spectrum, this efficiency was evaluated to be about 
72%. The reduced capture time window has an effi- 
ciency of 94%. The residual accidental backgrounds 
was thus reduced to < 0.05/day at the near sites, and 
< 0.01/day at the far site. 

To reduce the number of minimum ionizing muons 
in these data samples, we assumed that most of the 
9 Li and 8 He production was accompanied with neu- 
tron generation, and thus rejected AD tagged muon 
events with no follow-on neutron (defined as > 1.8 
MeV signal within a 10 - 200 /is window). The muon 
samples with and without reduction were both pre- 
pared for the 9 Li and 8 He background estimation. 
The data was sub-divided into six groups in visible 
muon energy (0.02-0.5, 0.5-1.5, 1.5-2.5, 2.5-3.5, 3.5- 
4.5, and > 4.5 GeV). Taking EH1 as an example, the 
corresponding muon rates in each energy bin were 
(10.0, 10.9, 0.23, 0.042, 0.016, 5.6e-3 Hz). Note that 
the maximum visible energy was around 5 GeV be- 
cause of the saturation of the PMTs. An example 
of a fit to the time-since-last muon distribution us- 
ing Eq. [8] for determining the number of 9 Li and 8 He 
events for > 4.5 GeV is shown in Fig. [T5] Though 
only four seconds are shown in the figure, the fit range 
was actually from 1 ms to 40 s. Fitting over such a 
large range helped to insure that i? M was accurate. 
Because of the 1000 /js had veto, the fitted i? M was 
slightly smaller than the directly measured value. 

« 
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Fig. 15. An example of fitting for 8 He/ 9 Li 
backgrounds. 

Instead of allowing the 9 Li to 9 Li plus 8 He ratio to 
float, we scanned it from to 1 in steps of 0.01. For 
each scan point, a maximum likelihood fit was done, 
where only N Li+He , N IBD , and i? M were allowed to 
float. Also, only the results with a global maximum 
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likelihood among scan points were regarded as best 
fit values. The global maximum likelihood confirmed 
that 9 Li was dominant in the 8 He/ 9 Li backgrounds. 
The binning effect was included in the uncertainty 
estimation by changing the bin width of the time- 
since-last muon distribution. 
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Fig. 16. The fitted 9 Li yield as a function of the 
visible energy of parent muons for three exper- 
imental halls. The open circles represent the 
fit with all muons included. Due to high muon 
rate, the fit is done only for > 2.5 GeV. 
The filled circles are the results obtaining by 
requiring a neutron following the muon as de- 
scribed in the text. In the bottom panel the 
prediction from the near site measurements is 
shown as solid line. 



The best-fit results are shown in Fig. [TBI Since the 
statistics were quite low in EH3, we also predicted the 
9 Li yield in EH3 from the EH1 and EH2 yields by as- 
suming that the 9 Li yield with the same visible energy 
at different sites were identical, as shown in the bot- 
tom panel in Fig. 1161 The measured values agreed 
with the prediction within statistics. Another check 
was done by predicting the EH3 9 Li yield assuming 
that it follows an E°- 74 power law, where is the 
simulated average muon energy (See Table [J), and 
normalizing to EH1 and EH2 measurements. Again 
the fitted EH3 9 Li yield agree with the prediction 
within statistics. By considering binning effects, dif- 
ferences between the results with and without muon 
reduction, and the difference between the predicted 
EH3 result and the measured result, we assigned a 
50% systematic uncertainty to the final result. 



5.4 (a, n) Backgrounds 

The 13 C(a, n) 18 background was determined by 
measuring alpha-decay rates in situ and then using 
the MC to calculate the neutron yield. We identified 
four sources of alpha decays, the 238 U, 232 Th, 227 Ac 
decay chains and 210 Po. The decay chains are f3-a 
cascades with half lives of 164.3 /is, 0.3 /is, and 1.781 
ms, respectively. Fig. [T7] displays the correlation of 
the prompt-delayed energy distributions for various 
time intervals corresponding to these cascade decays: 
1-3 /xs at upper left (group A are 212 Bi- 212 Po decays 
from the 232 Th decay chain), 10-160 [is at upper right 
(group B are IBD events where the neutron captures 
on hydrogen. Group C are 214 Bi- 214 Po decays from 
the 238 U decay chain, and group D are 219 Rn- 215 Po 
decays from the 227 Ac decay chain). In the 1-2 ms re- 
gion at lower left, only group D and some accidental 
coincidence events remain. 




Delayed Energy (MeV) 



Delayed Energy (MeV) 



Fi. 



g. 17. Correlations of prompt and delayed en- 
ergy for cascade decay chains of contaminants 
within the ADs. At upper left are events with 
a time correlation between 1 and 3 [is, 10 to 
160 fis is at upper right, 1-2 ms is at lower 
left, with the combined distributions at lower 
right. 



210 Po was produced by the decay of 222 Rn. Its 
5.3 MeV alpha produced 0.5 MeV of visible energy 
in an AD. The spatial distribution suggests that the 
210 Po background was due in part to an accumulation 
on the wall of the inner acrylic vessel. 

Geant4 was used to model the energy deposition 
process. Based on the (a, n) cross sections archived 
in JENDL [321 . the neutron yield as a function of 
a energy was calculated and summed. Finally, with 
the in-situ measured alpha-decay rates and MC de- 
termined neutron yields, the 13 C(a, n) 16 rate was 
calculated, as listed in Table [5] The uncertainties 
come from the selection efficiencies of the Bi-Po and 
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Rn-Po chain measurements, the possible deviation 
from equilibrium of the 238 U, 232 Th, and 227 Ac de- 
cay chains, the fitting to determine the 210 Po activ- 
ity, and the simulation of (a, n) reactions. During the 
Gd-LS synthesis, 238 U, 232 Th, and Ra were removed 
by radio-purification. They may contribute *~30% of 
the alphas of the whole chain. Thus a 30% uncer- 
tainty was assigned for the possible deviation from 
equilibrium, which was the largest component in the 
uncertainties. A 10% uncertainty was assigned to the 
neutron yield by comparing the MC simulation with 
an analytical calculation. Together with the other 
two components, ^50% uncertainties were estimated 
for the (a, n) backgrounds, slightly different for each 
AD due to different alpha components in them. 

5.5 Correlated Backgrounds from Am-C 
Source 

During data taking, the Am-C sources sat inside 
the ACUs on top of each AD. Neutrons emitted from 
these sources would occasionally mimic IBD events 
by scattering inelastically with nuclei in the shield- 
ing material (emitting gamma rays) before being cap- 
tured on a metal nuclei, such as Fe, Cr, Mn or Ni 
(releasing more gamma rays). It was possible for the 
gamma-rays from both processes to enter the scin- 
tillating region and satisfy the IBD selection require- 
ments. Fig. llSl shows the energy spectrum in the three 
ADs at the far site of these delayed candidates from 
the Am-C sources. The rate in MC was normalized 
to data. There is good agreement between the data 
and MC. 




Energy (MeV) 



Fig. 18. Energy spectrum for events near the 
top of the three ADs in the far hall show three 
peaks consistent with neutron capture on 56 Fe 
and 58 Ni/ 54 Fe/ 53 Cr. 

Fig. [19] shows an asymmetry of delayed- type 
events along the z axis as was seen by ADs in the far 



hall. We estimated the delayed-type events from the 
Am-C sources by subtracting the number of delayed- 
type singles in the Z < region from the Z > re- 
gion. The Am-C correlated background rate was esti- 
mated by MC simulation normalized with the Am-C 
delayed-type event rate obtained from the data, 



J-^n—HheMC 

where N corr - M c and N n -uh e uc are the number of 
correlated background and number of delayed-type 
events in MC respectively, and R n -uke data is the Am- 
C delayed-type event rate from data. Even though 
the agreement in shape between data and MC is ex- 
cellent for Am-C delayed-type events, we assigned 
100% uncertainty to the estimated background due 
to the Am-C sources to account for any potential un- 
certainty in the neutron scattering/capture cross sec- 
tions used in the simulation. 




-3-2-10 1 2 3 

Z(m) 

Fig. 19. Z distribution of delayed-type events. 
The excess in the top half of the ADs (Z > 0) 
comes from the Am-C sources in the ACUs. 

6 Side-by-Side Comparison in EH1 

Relative uncertainties were studied with data by 
comparing two side-by-side antineutrino detectors. A 
detailed comparison using three months of data from 
the ADs in EH1 has been presented elsewhere jl5j |. 
An updated comparison of the prompt energy spec- 
tra of IBD events for the ADs in EH1 using 231 days 
of data (Sep. 23, 2011 to May 11, 2012) is shown 
in Fig. [20] after correcting for efficiencies and sub- 
tracting background. A bin- by-bin ratio of the AD1 
and AD2 spectra is also shown. The ratio of the to- 
tal IBD rates in AD1 and AD2 was measured to be 
0.987±0.004(stat.)±0.003(syst.), consistent with the 
expected ratio of 0.982. The deviation of the ratio 
from unity was primarily due to differences in the 
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baselines of the two ADs with a slight dependence 
on the individual reactor on/off status. It was shown 
that AD2 has a 0.3% lower energy response than AD1 
for uniformly distributed events, resulting in a slight 
tilt to the distribution shown in the bottom panel of 
Fig. [20] The distribution of the data points denoted 
by open circles was created by scaling the AD2 en- 
ergy by 0.3%. The bin-by-bin ratio with scaled AD2 
energy agrees well with a flat distribution. 
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Prompt Energy (MeV) 
Fig. 20. The energy spectra for the prompt 
signal of IBD events in AD1 and AD2 are 
shown in the top panel, along with the bin-by- 
bin ratio in the bottom panel (solid circles). In 
the bottom panel, the dashed line represents 
the ratio of the total rates for the two ADs, 
and the open circles show the ratio with the 
AD2 energy scaled by +0.3%. 

7 Reactor Antineutrino Flux 

Reactor antineutrinos result primarily from the 
beta decay of the fission products of four main iso- 



topes, 



5 U, 239 Pu, 238 U, and 241 Pu. The D e flux of 



each reactor (S(E)) was predicted from the simulated 
fission rate (Fj) a nd the antineutrino spectrum per 
fission (Si) 33rl38( of each isotope [39(, 

S(E) = Y,F i S i (E) (10) 

i 

where i sums over the four isotopes. The fission rate 
was determined from the fission fraction the en- 
ergy released per fission E il and the measured ther- 
mal power (W t h)) 



where both i and k are indices over the four isotopes. 

The thermal power data were provided by the 
power plant. The uncertainties were dominated by 
the flow rate measurements of feedwater through 
three parallel cooling loops in each core 3914411]. The 
correlations between the flow meters were not clearly 
known. We conservatively assumed that they were 
correlated for a given core but uncorrelated between 
cores, giving a maximal uncertainty for the experi- 
ment. The assigned core-to-core uncorrelated uncer- 
tainty for thermal power was 0.5%. 

A simulation of the reactor cores using commer 
cial software (SCIENCE " 



42, 4 



provided the fission 
fraction as a function of burn-up. One example of 
fuel evolution is shown in Fig. [21] The fission frac- 
tion carries a 5% uncertainty set by the validation of 
the simulation software. A complementary core simu- 
lation package was developed based on DRAGON [44 [ 
as a cross check and for systematic studies. The code 
was validated with the Takahama-3 benchmark [45| 
and agrees with the fission fraction provided by the 
power plant to 3%. Correlations among the four iso- 
topes were studied using the DRAGON-based simula- 
tion package, and agrees well with the data collected 
in Ref. [46(. Given the constraints of the thermal 
power and correlations, the uncertainties of the fis- 
sion fraction simulation translated into a 0.6% core- 
to-core uncorrelated uncertainty in the neutrino flux. 
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Fig. 21. Fission fractions of reactor fuel iso- 
topes as a function of burn-up from a simula- 
tion of reactor core Dl. Other isotopes con- 
tribute less than 0.3% in total. 

The antineutrino spectra per fission is a corre- 
lated uncertainty that cancels out for a relative mea- 
surement. The reaction cross section for isotope i 
was defined as rjj = J S i (E JJ )a(E v )dE l/ , where Si(E v ) 
is the antineutrino spectra per fission and a(E v ) is 
the IBD cross section. We initially took the reac- 
tion cross section from Ref. but substituted the 
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IBD cross section with that in Ref. [48|. The en- 
ergy released per fission and its uncertainties were 
taken from Ref. [i^ . Non-equilibrium corrections for 



long-lived isotopes were applied following Ref. 37 1. 
Contributions from spent fuel 50|, l5l| (~0.3%) were 
included as an uncertainty. 

Table 6. Reactor-related uncertainties. 





Correlated 




Uncorrelated 




Energy /fission 
IBD reaction/fission 


0.2% 

3% 


Power 
Fission fraction 
Spent fuel 


0.5% 
0.6% 
0.3% 


Combined 


3% 


Combined 


0.8% 



The 3D spatial distribution of the isotopes within 
a core was also provided by the power plant. Sim- 
ulation indicated that the spatial distribution has a 
negligible effect. As such, the reactor core was taken 
as a point source. The 3D core simulation with the in- 
put of the monthly in-core neutron flux measurement 
showed that the fission gravity center moves less than 
1 cm on the horizontal plane and several cm vertically 
as the fuel burned. The resulting baseline variation 
can be ignored. 
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Fig. 22. The measured daily average IBD rates 
per AD in the three experimental halls are 
shown as a function of time along with pre- 
dictions based on reactor flux analyses and 
detector simulation. The reactor status is in- 
dicated. 

Fig. [22] presents the background-subtracted and 
efficiency-corrected IBD rates in the three experi- 
mental halls. The predicted IBD rates from reactor 
flux calculation and detector simulation are shown for 
comparison. The dashed lines have been corrected 
with the best-fit normalization parameter e in Eq. [13] 
to reduce the biases from the absolute reactor flux 



uncertainty and the absolute detector efficiency un- 
certainty. These predictions are systematically higher 
than the data points due to the oscillation effect in the 
data at the near and far sites. Predictions accounting 
for oscillation effects are also shown. 

8 Results 

The V e rate in the far hall was predicted with a 
weighted combination of the two near hall measure- 
ments assuming no oscillation. A ratio of the mea- 
sured to expected rate is defined as 



R-- 



Mj_ 



M f 



aM a + fiM b 



(12) 



where N j and Mf are the predicted and measured 
rates in the far hall (sum of AD 4-6), M a and M b 
are the measured, background-subtracted IBD rates 
in EH1 (sum of AD 1-2) and EH2 (AD3), respec- 
tively. The weights a and /3 are not unique since 
we approximate the contributions of the six reactors 
with the two observables M a and M b . All valid phys- 
ical models to determine the weights should satisfy 
the normalization requirement, i.e. the combination 
of the two near hall predictions should be equal to 
the direct prediction of the rate in the far hall in 
terms of antineutrinos emitted by the reactors. In 
this analysis we also required a maximum cancella- 
tion of reactor uncertainties, ignoring the statistical 
and detector-related systematic uncertainties. The 
values for a and (3 were dominated by the baselines, 
and only slightly dependant on the integrated flux of 
each core. For the analyzed data set, a = 0.0444 and 
(3 — 0.2991. The residual reactor-related uncertainty 
in R was 5% of the uncorrelated uncertainty of a sin- 
gle core. The ratio observed at the far hall was: 

R = 0.944 ± 0.007(stat.) ± 0.003(syst.) , 

where the statistical (systematic) uncertainties were 
obtained by propagating statistical (uncorrelated sys- 
tematic) uncertainties in the measured IBD counts in 
the three halls. 

The value of sin 2 2#i 3 was determined with a \ 2 
constructed with pull terms accounting for the corre- 
lation of the systematic errors , 



E 



[M d -T d (l+e + J2 r ^a r +e d )+7 ld \ 2 



M d + B d 



(13) 
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where M d are the measured IBD events of the d-th 
AD with its backgrounds subtracted, B d is the corre- 
sponding background, T d is the prediction from an- 
tineutrino flux, including MC corrections and neu- 
trino oscillations, u>^ is the fraction of IBD contribu- 
tion of the r-th reactor to the d-th AD determined 
by the baselines and antineutrino fluxes. The un- 
corrected reactor uncertainty is o r (0.8%), as shown 
in Table [S] The parameter a A (0.2%) is the uncor- 
rected detection uncertainty, listed in Table |4j The 
parameter a B is the quadratic sum of the background 
uncertainties listed in Table [5j The corresponding 
pull parameters are (a r ,e dl r] d ). The detector- and 
reactor-related correlated uncertainties were not in- 
cluded in the analysis. The absolute normalization e 
was determined from the fit to the data. 

The survival probability used in the \ 2 was 



1 - sin 2 20 13 sin 2 : (1.267 Am 2 31 L/E) 
cos 4 13 sin 2 20 12 sin 2 (1 .267 Am 2 21 L / E) , 



l 31 



' 2 ,sin 2 20 12 = 0.861™, 



and Am 2 j = 7.59±%£ xlO" 5 eV 2 



53j . The uncertainty 



in Am,! [lj| had negligible effect and thus was not 
included in the fit. 
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Fig. 23. Ratio of measured versus expected 
signals in each detector, assuming no oscilla- 
tion. The error bar is the uncorrelated un- 
certainty of each AD, including statistical, 
detector-related, and background-related un- 
certainties. The expected signal has been 
corrected with the best-fit normalization pa- 
rameter. Reactor and survey data were used 
to compute the flux-weighted average base- 
lines. The oscillation survival probability at 
the best-fit value is given by the smooth curve. 
The AD4 and AD6 data points were displaced 
by -30 and +30 m for visual clarity. The \ 2 
value versus sin 2 2f?i3 is shown in the inset. 



The best-fit value is 

sin 2 20 1Z = 0.089 ± 0.010(stat.) ± 0.005(syst.) 

with a x 2 /NDF of 3.4/4. All best estimates of pull pa- 
rameters are within its one standard deviation based 
on the corresponding systematic uncertainties. The 
no-oscillation hypothesis is excluded at 7.7 standard 
deviations. Fig. [53] shows the number of IBD can- 
didates in each detector after correction for relative 
efficiency and background, relative to those expected 
assuming no oscillation. A ~1.5% oscillation effect 
appears in the near halls, largely due to oscillation of 
the antineutrinos from the reactor cores in the far- 
ther cluster. The oscillation survival probability at 
the best-fit values is given by the smooth curve. The 
X 2 value versus sin 2 2f? 13 is shown in the inset. 

The observed T7 e spectrum in the far hall was com- 
pared to a prediction based on the near hall measure- 
ments aM a + j3M b in Fig. [Ml The distortion of the 
spectra is consistent with that expected due to oscilla- 
tions at the best-fit 13 obtained from the rate-based 
analysis. 
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Fig. 24. Top: Measured prompt energy spec- 
trum of the far hall (sum of three ADs) com- 
pared with the no-oscillation prediction based 
on the measurements of the two near halls. 
Spectra were background subtracted. Uncer- 
tainties are statistical only. Bottom: The ra- 
tio of measured and predicted no-oscillation 
spectra. The solid curve is the expected ra- 
tio with oscillations, calculated as a function 



of neutrino energy assuming sin 26*i 



0.089 



obtained from the rate-based analysis. The 
dashed line is the no-oscillation prediction. 
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9 Conclusions 

We have updated the measurement of the neu- 
trino mixing angle 9 13 with a 116.8 kton-GW th -day 
livetime exposure at the far hall. A total of 138,835, 
66,473, and 28,909 electron antineutrino candidates 
were detected in the Daya Bay near hall, the Ling 
Ao near hall, and the far hall, respectively. Com- 
pared with the prediction based on the near-hall mea- 
surements, a deficit of 5.6% was observed in the far 
hall. The rate-based analysis has yielded sin 2 2(9 13 = 
0.089 ± 0.010(stat.) ± 0.005(syst.). This is the most 
precise measurement of sin 2 2# 13 to date with a pre- 
cision of 12.6%, and supersedes our previous mea- 
surement . We anticipate additional improvements 
following the installation of two additional ADs in 
advance of an extended data run. 
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